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We demonstrate that the Bell test cannot be realized at finite temperatures in the vast majority of electronic 
setups proposed previously for quantum entanglement generation. This fundamental difficulty is shown to 
originate in a finite probability of quasiparticle emission from Fermi-sea detectors. In order to overcome the 
feedback problem we suggest a detection strategy, which takes advantage of a resonant coupling to the quasi- 
particle drains. Unlike other proposals the designed Bell test provides a possibility to determine the critical 
temperature for entanglement production in the solid state. 
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It is well-known that, unlike photons, quasiparticles in the 
Fermi sea injected from reservoirs, which are kept at thermal 
equilibrium, can be entangled by just a tunnel barrier. 1 This al- 
lows for particularly simple proposals for quantum quasipar- 
ticle entanglement, which do not involve interactionS i 1 ' 2 ' 3 ' 4 ' 5 
Theoretical results for the entanglement production in differ- 
ent electronic setups have been summarized in Refs.|gand|7l 
while yet no experimental evidence of the quasiparticle entan- 
glement in the Fermi sea has become available. 

The quantum entanglement of two particles with respect to 
a spin-like degree of freedom can be accessed experimentally 
by measuring the spin correlator 



C(a,b) = ((a-a) 1 ®(h-(T) 2 ), 



(1) 



where er = (a x , a y ,a z ) is the vector of Pauli matrices. The 
spin projection of the particles in the detectors 1 and 2 is mea- 
sured with respect to the unit vectors a and b, correspond- 
ingly. If the correlation between the particles is of a classical 
origin the following Bell inequality holds^ 

B = \C(a, b) + C(a', b) + C(a, b') - C(a', b')| < 2, (2) 

for arbitrary choice of the unit vectors a, b, a', b'. The vi- 
olation of the inequality (|2]i is, therefore, sufficient but not 
necessary condition for quantum entanglement. 

In solid-state electronics we deal with elementary excita- 
tions in the Fermi gas, which are referred to as quasiparti- 
cles. Even though the pairwise quasiparticle entanglement^^ 
is believed to be generated in many deviceS j n i 12 i 13 i 14 its ex- 
perimental observation is obscured by the nature of electronic 
detectors. Those, unlike the photodetectors in optical setups, 
contain a number of quasiparticles in the ground state, which 
fill up available quantum levels below the Fermi energy. If 
a part of the device is at finite temperature the electron and 
hole excitations are spontaneously created near the Fermi sur- 
face resulting in a finite probability for a Fermi-sea detector 
to emit. Such processes are harmful for any sensible Bell test. 

The problem of quasiparticle entanglement detection has 
been put forward in Refs.ll5landll61 where the possibility to 
construct a Bell-type inequality with current cross-correlators 
is discussed. It has been suggested to take advantage of the 



generalized spin-correlator 

((/V 1T 



C M (a,b) 



N U ) (7V 2T - N 2l )) 



((JV 1T + N n ) (N 21 + N 2l )) ' 



(3) 



where N ncr is a number of particles with a spin projection a 
registered by the detector n. (In solid state the role of spin 
can be played by other quantum degrees of freedom such as 
orbital momentum or isospin). Similarly to Eq. (Q~|i the spin 
projection in Eq. (0 is measured with respect to the direction 
a in the first detector and b in the second one. Both definitions 
(Q]l and ® are equivalent in the original Bell setup, if no more 
than two particles are received within the detection time and 
the detectors do not emit particles. In electronic circuits the 
number of quasiparticles N na is given by the time integral of a 
current I na flowing to the corresponding Fermi-sea reservoir 



N„ 



oc 



tdct 

dtl na (t), 



(4) 



which is not restricted. For large detection times td e t one typ- 
ically observes |7Y n -f — N n ± \ <C |JV n | + N n i\, hence the Bell 
inequality (|2]) cannot be violated and the corresponding mea- 
surement is useless for entanglement detection. The difficulty 
has been discussed in Ref.[l6| for zero temperature. 

An essential problem occurs in the opposite limit tdct — ► 0, 
because N na defined by Eq. © can take on negative values. 
This leads to fluctuations with |JV n f — N n ±\ > \N n -\ + iV n j.|> 
which are explicitly forbidden in the Bell test. This situation is 
realized at finite temperatures. Then, the violation of Eq. (0 
has no relation to the entanglement detection and the corre- 
sponding measurement is not of a Bell type. 

Thus, the violation of the inequality (|2]i with the corre- 
lator C substituted by C M does not provide a conclusive 
evidence for quantum entanglement generation at any fi- 
nite temperature. This difficulty clearly applies to the de- 
tection of electron-hole entanglement 713 produced by tun- 
neling events or by time-dependent gating. But even in 
more sophisticated setups where zero-temperature detectors 
and finite-temperature sources are represented by different 
metallic leads (in close resemblance to the original Bell pro- 
posal) the Bell test based on Eq. ([3]) is flawed. Exam- 
ples include three-terminal fork geometries 4 and four-terminal 
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beam-splitter geometries with grounded detectors. We focus 
on the latter (see Fig. Q]) due to a number of previously pro- 
posed realizations , 'i 3 ^ 5 ' 11 ' 17 ' 18 which are mostly based on the 
directed transport along quantum-Hall edge channels. Minor 
modifications, such as lowering chemical potential in one of 
the detectors or increasing detection time, can suppress the 
probability of detector emission but lead, instead, to useless 
measurement with | N n ^ — N n i | <C | N n -[ + N n i | . The generic 
situation is illustrated in Fig.[2]for the case of electronic beam 
splitter. 

For quantum particles the spin correlator from Eq. ([3]) is 
expressed through the expectation value 



(N la N 2(7 ,) cx K a 



(5) 



l ^aa' — <-det 



{ha)(h,')+ I —V aa ,{u)F{ut Aet /2) 



where N and / are regarded as operators. We introduce the 
function T(x) = (sin x) 2 jx 1 and the frequency-dependent 
cross-correlator 



= / dte iut (6h a (t)I 2(T ,(0)), 



(6) 



with SI na (t) = I n<T (t) — (I na )- In Figs. 12131 the correlator 
C M defined by experimentally measurable quantities ( 15l6t is 
compared with the exact result of the density matrix analysis 
of the final stated 

We restrict ourselves to an important class of systems which 
do not involve spin-dependent scattering, because the chances 
to generate quantum entanglement with respect to the spin de- 
gree of freedom are obviously maximized in such setups. The 
values of V aa ' in Eq. © are related to the cross-correlator V 
of the corresponding spin-independent problem as 



7>Tt =Vu = |(l + ab)7>, 
T U =V U = §(l-ab)P. 



(7) 
(8) 



This symmetry holds even for interacting electronic systems 
provided the absence of spin dephasing. It follows from 
Eqs. ( 1718b that a neglection of the mean currents 5 - 7 in Eq. (J5]l 
is equivalent to C M (a, b) = ab, hence the inequality (O is 
violated with £> max = 2\/2 irrespective of voltages, tempera- 
ture or other setup characteristics. Clearly, such violation has 
nothing to do with pairwise quantum entanglement. We will 
see that the problem persists even if the exact expression for 
K aa i is used. 

In the absence of spin-dependent scattering the mean cur- 
rents measured by the detectors do not depend on the direc- 
tions a and b, (I na ) — (In)- From Eqs. (I5I6I7I8I ) we obtain 



C M (a,b) 



7 



■ ab, 



2 + 7 

with a parameter 7 given by the ratio 

/ du> V{u)T{ujt dct /2) 



1 = 



Mh){h) 



(9) 



(10) 



where both the cross-correlator V(u) and the product of the 
mean currents have to be calculated for the corresponding 
spin-independent problem. 
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FIG. 1 : A generic beam splitter for entanglement production in the 
solid state. The voltage bias applied between the sources SI and 
S2 generates an entangled outgoing state at the detectors Dl, D2 
provided the temperature in the sources T is smaller than a critical 
temperature T c . 



An example of such a calculation can be performed within 
the Landauer-Biittiker scattering approach, 19 which is valid 
as far as inelastic processes in between the reservoirs can be 
disregarded. Within the scattering approach the mean current 
to the reservoir a is given by 



(Ia) = 



JdEj2 (<W - \S a p{E)\ 2 ) ME), (11) 



where f a (E) = (1 + exp [(E - eV a ) / k^Ta])' 1 is the Fermi 
distribution function, which depends on the temperature of the 
corresponding reservoir T a and the voltage bias V a applied. 
The frequency-dependent correlator © of the currents flow- 
ing to the reservoirs a and a' reads^ 

e 2 f . , 

V aa '(uj) = — / dE y j M aa i£p, (E, huj)F j3 p l (E, hu>), 
J pp 

F^,(E,Q) = U(E)fp(E + + fp{E)fp,{E + Q), 
M aa ,,^{E,n) = (5 a0 5 a0l - S*p(E)S af3 ,(E + «)) 

X {SafpScp - S* a , p/ (E + n)Sa>p{E)) , (12) 



f(E) ee 1 - f(E). 



(13) 



Let us consider a generic beam splitter with no spin- 
dependent scattering depicted schematically in Fig. Q] Such 
a setup is characterized by an energy-independent 5-matrix 



S 



s' 
s 



(14) 



where 2x2 unitary matrices s and s' describe the transport 
from sources to detectors and from detectors to sources, cor- 
respondingly. We parameterize 



o}4> 






,i<j>' 



\/l - T iy 











(15) 



where r € [0, 1] is the beam-splitter transparency and the 
spin index is omitted. Following the majority of propos- 
als both detectors and the second source are grounded, i.e. 
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FIG. 2: The spin cross-correlator C obtained from the density matrix 
of the final scattering state (solid lines; cf. Eg. 120b. and its general- 
ization C M , evaluated numerically from Eqs. J9ll5b for different val- 
ues of the detection time eVt dct /h = 0.01(»), 0.1(a), 1(b), 5(») 
(dashed lines; see Eqs. 121 1221 . 



Vb = Vm = Vd2 = 0, V S2 = 0, while V S i = V is the 
voltage applied between the sources. 

At zero temperature the beam splitter acts 7 as a source of 
spin-entangled Bell pairs 

|*b> = ^=1 Tiia - iiTa) (16) 

where the index n — 1,2 refers to the detector number. 
Such an entanglement generation is due to the Pauli princi- 
ple, which guarantees that a filled state with E G (0, eV) in 
the first source contains exactly two quasiparticles with the 
opposite spins. 

The Bell pairs can be accessed at zero temperature by per- 
forming a time coincidence detection. For finite temperature 
T in the sources the density matrix projection, which corre- 
sponds to a single particle in each detector, is derived in the 
AppendixfAf 

where I4 is the unit matrix in the two-particle Hilbert space 
and £ is an energy-independent weight factor 

{ = i^ii^f^l ^ (18) 

r(l-T)(fsi-fs2) 2 + 2fsifsifs2f S 2 

where f$ n is the Fermi distribution function in n-th source. 
The result ( fT7b describes the mixed Werner stated which 
is entangled as far as £ > 1/3 according to the Wootters 
formula. 21 In the present case this condition is equivalent to 
T < T c with the critical temperature T c determined by the 
equation 7 - 

r(l - r) sinh 2 (eV r /2fc B T c ) = 1/4. (19) 
From Eqs. ( 111171 ) one obtains the exact spin correlator 

C(a,b) = -£ab, (20) 



which is plotted in Fig. [2] with the solid line for different val- 
ues of the transparency parameter. The corresponding Bell 
inequality Q can be violated for £ > 1/V2, which is, in- 
deed, a sufficient condition for the entanglement. Whether or 
not such a Bell test can be performed by measuring current 
cross-correlator (0 is, however, an open question. 

In order to answer this question we substitute the expression 
(TBI for the 5-matrix to Eqs. dl 1I121 >. where the summation 
runs over the index a = {SI, S2, Dl, D2}. The correlator 
C M is, then, obtained from Eqs. ( 191101 ) with h = I m , I 2 = 
I D2 , andV = V D \,D2- 

For tdet 3> mm{h/eV, h/k#T} we obtain 

i.e. the corresponding measurement is useless for an entangle- 
ment detection. Indeed, such a long-time measurement is not 
projective, therefore it does not single out the state with one 
quasiparticle in each detector^ 

In the opposite limit we, however, find 

7 = -1, <dct < mm{h/eV, h/k B T}, (22) 

hence the inequality (O is violated for any temperature of the 
source. Thus, according to the density matrix analysis dl7|20K 
the corresponding measurement is not of a Bell type. Both 
results (|2TT > and d22l formally hold for any temperature of the 
detectors. 

The transition from non Bell-type measurement to the use- 
less measurement with the increase of idot is illustrated in 
Fig. [2] The Bell parameter defined with the correlator C M 
does not depend on the beam-splitter transparency r and can 
easily exceed 2 even in the absence of any entanglement. 

The result of Eq. d22l) is equivalent to 

(lDi(t)W*)}=0. (23) 

At T = the currents Io n (i) are sign-definite, hence Eq. (l23l 
is exact for every single time-coincidence measurement in 
agreement with the prediction of the density matrix approach. 
For rising temperatures T > the correlation (l23l holds only 
on average and is not sensitive to vanishing quantum entangle- 
ment in the final state of the beam splitter (fTTT i. Consequently, 
the inequality (O with C substituted by C M can be violated for 
arbitrarily high temperatures. The absence of critical temper- 
ature indicates once again 23 that such a violation has nothing 
to do with the entanglement detection. Instead, the decay of 
C M (a, b) with the temperature in Fig. [2] (dashed lines) is de- 
termined by the detection time t^et- 

Thus, the measurement of C M (a, b) cannot be used for the 
entanglement test in the beam-splitter setup and the value of 
T c cannot be inferred from such a measurement as the matter 
of principle. 

We propose a way to rescue the Bell measurement by cou- 
pling detectors via the energy filters, which are described by 
energy dependent scattering amplitudes: r n ,r' n , t n , t' n , where 
n = 1, 2 is the number of the detector. The use of energy- 
filters in the context of Bell measurement at zero temperature 
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FIG. 3: The case of resonant detector coupling. The short-dashed 
line shows C M from Eq. l !26t . while the long-dashed lines are nu- 
merical results for Breit-Wigner resonances J24b with finite width 
F = O.OleV, detector voltage Vd ~ —V, and different detection 
times Ttdct/h — 0.01(»), 0.1(A). The measurement is useless for 
idet > 0.1/i/r. The solid line shows the correlator C from Eq. J20t . 

has been discussed in Ref. 24. Let us illustrate our results for 
the case of identical filters with the Breit-Wigner form of the 
transmission amplitude 

t„(E) = e lS "(Y/2)(E -E - ^^/2)- 1 . (24) 

The 5-matrix of the full setup including the filters is given by 

( s'r{E)s s't'(E)\ 
S{E > ~ { t(E) s r'(E) ) (25) 

where t = diag(ti, t2), r — diag(ri,r2), etc. The condi- 
tion for time-coincidence detection now reads tdct <C h/Y. 
The currents Idu (t) can be made sign-definite by applying an 
additional voltage bias Vd, as shown in Fig. Q] The current 
fluctuations due to temperature are not harmful for the Bell 
test as far as \bVd\ 3> k B T, which is the only restriction on 
the value of Vd- In this case there is no requirement for an 
additional cooling of the detectors, so that a whole setup can 
be kept in temperature equilibrium. Moreover, for Y <C eV 
the dependence on t^ct vanishes, meaning that C M (a, b) can 
be obtained experimentally from zero-frequency noise mea- 
surements. The feasibility of such a Bell test is illustrated in 
Fig. [3]for realistic values of the parameters. 
For T -> we obtain from Eqs. ( I9ll0ll2l24l25l l 

c m = _ r(l - r)(f sl - f S2 ) 2 ah 

r(l-r)(fsi-fs2) 2 +2fsifs2 E=Eo ' 

The result is plotted with the short-dashed line in Fig. [3] It is 
evident from the comparison with Eqs. (118120b that the pro- 
posed measurement is always of the Bell type. The corre- 
lator ( f26b tends to the exact one (l20l for Eq 3> eV. The 
setup efficiency is, however, exponentially low in this limit. 
The numerical results in the case of finite resonance width 
T = O.OleV^ are plotted in Fig.|3]with the dashed lines. The 
test provides the lower estimate for the critical temperature. 



In conclusion we point out the fundamental restrictions for 
the Bell test in electronic setups due to the quasiparticle emis- 
sion from Fermi-sea detectors. We propose a way to rescue 
the Bell measurement by a resonant coupling to the detectors. 
We show that the lower estimate of the critical temperature 
for entanglement production can be experimentally obtained 
in the proposed setup. 

This research was supported by the DFG Priority Pro- 
gramm 1285. The discussions with C. W. J. Beenakker, 
W. Belzig, and Yu. V. Nazarov are gratefully acknowledged. 

APPENDIX A: DENSITY MATRIX PROJECTION 

Following Ref. @ we review the derivation of the density 
matrix projection ( [171 ) for final scattering state in the case of 
the setup depicted in Fig.Q] The density matrix of the incom- 
ing state is given by 

P in = LI (fsn(E)\0)(0\+fsn(E)ai aE \0)(0\a naE ), 

n,E,a 

(Al) 

where fs n is the Fermi distribution function in the source S n , 
fsn = 1 — fsn, and a naE is the fermion annihilation operator 
for an incoming scattering state at the channel n and energy E. 
The annihilation operators for the outgoing scattering states, 
b naE , are obtained from the relation, 

ma h, i 

(A2) 

m 

where s nm are the components of a unitary scattering matrix. 
Thus, the density matrix of the final state is 

p out = JJ ^ Sn{E) |0)(Q| + fsn[E) c t^ | )(0| C „ CTB } , 
n . E . a 

(A3) 

where 

cLE=Y, b l°ESmn{E). (A4) 

m 

In order to quantify the two-particle entanglement for the 
partition H.D1 <8> "Hdi of Hilbert space with respect to the 
detectors, the state p ont has to be projected onto the sec- 
tors M Ei n 1 ,e 2 n 2 of the Fock space with the energies E\, E 2 
and particle numbers N\ , N 2 in the corresponding detectors 
Di, D 2 - The density matrix No of the projection factor- 
izes into a product state in all sectors except for the sector 
Nei,ei with E\ = E 2 = E and ^1=^2 = 1. Projection 
onto this sector is found from Eq. ( IA3b as 

W\\p\\ — fsifsifS2fs2 I4 

+ 2r(l-r)(/ sl - fs2) 2 \^B){^ B \, (A5) 

where I4 is the unit matrix in the two-particle Hilbert space, 
\^b) is the Bell state ( [To*] ), and the weight factor wu is deter- 
mined from the condition Tr pn — 1 as 

wu = ifsifsifs2fs2 + 2t(1 - r)(/ S i - f S 2) 2 - (A6) 
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From Eqs. ( IA5IA6I ) we obtain Eqs. ( 117118b . By substituting 

fsn = (1 + exp[(E - eV n )/k B T]) with V 1 =V,V 2 = in 
Eq. (fT8l we can further simplify the parameter £ as 



£CO = i - 



1 + 2r(l - r) sinh 2 



eV 
2k B T 



(A7) 



The critical temperature T c is determined from the equation 
£(T C ) = 1/3, which is equivalent to Eq. dT9b . 



APPENDIX B: EVALUATION OF THE CORRELATOR C M 

1. Plain beam splitter 

We evaluate the generalized spin-correlator C M (a, b) given 
by Eqs. (19110b in the framework of the scattering approach. 
By substituting the scattering matrix ( TT4l Q3), into Eq. (fTTT i 
we calculate the mean currents, which are measured in the 
detectors Di,D 2 , as 



(Im) 



(1 - r) eV, (I D2 ) 



The cross-correlator (fT2l > is found as 



2fiw coth 



■reV. 



(Bl) 



C=±i 



2k B T 



(B2) 



The parameter 7 given by Eq. ([Tol l can be calculated analyti- 
cally in two opposite limits: 

(i) For large detection times, i<jet 3* min{/i/eV, /i/fceT}, 
one can replace £det-^"(wtdet/2) with 27T(5(w), hence 
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PC1,J2(0) 
*det (Idi)(Id2) 



(B3) 



This leads to the result (f2Tb . 

(ii) For short detection times tdct <C min{/i/ey, h/k B T} 
one can approximate !F(ojtd e t/2) ~ 1 in the relevant fre- 
quency range \ fhu\ < eV . In this limit the integral in Eq. ( TTOb 
does not depend on temperature 

j duVmM") = - 27r iff) r (! - T )> ( B4 > 



which leads to the simple result d22"i i. 



2. Beam splitter with energy filters 

We repeat the calculation in a more general case of an 
energy-dependent scattering matrix ( f25l l. From Eq. ( fTTT i we 
obtain the mean currents 



{I 



Dl 



l J dE\ tl (E)f 



dE\t 2 (E)f 



x lf D (E) - (1 - T)f sl (E) - Tf S2 (E)} , (B5) 

e 

x - rfsi(E) - (1 - r)/ S2 (£)] , (B6) 

where /d(-E) = (1 + exp [(E - eV D )/k B T}) is the Fermi 
distribution function in the detectors. From Eq. ( TT2b we cal- 
culate the cross-correlator 



ffll,D2(w) 



e 

2h 



r(l-r) | 



x i* (E) ti (E + fuJ) t 2 (E + hw) t 2 (E) 
x [F S i, S i(E,huj) + F S2 ^ S2 (E, fwj) 

-FsiME, M ~ F S2 ,si(E, M] , (B7) 

where the function F Q( 3 is defined in Eq. ( TT2b . These expres- 
sions allow for the numerical evaluation of C M (a, b) for arbi- 
trary energy-dependent scattering matrix $25[ . 

In the case of sharp resonances, such as those of the Breit- 
Wigner form d24l i with T — > 0, we have 



limT-^E)] 2 = -S(E - E ), 
and obtain from Eqs. ( 1B5IB6IB7IB81 ) 



(B8) 



7 = 



-r(l - r) (/ S1 - / S2 )' 



[/c- (1 - r)fsi- rf S2 ][fD- rfsi - (1 - r)f S2 ] ' 

(B9) 

where the Fermi functions are evaluated at the position of the 
resonance E = Eq > 0. For |eVb| 3> k B T the value Jd(Eo) 
is exponentially small, hence Eq. d26| i is justified. In general 
the setup is functional provided the detector voltage Vd is suf- 
ficiently large to ensure that fp is much smaller than both fsi 
an d fs2 within the energy window of the filters. 
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